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Abstract. In the past decade, most effort in the study of supermassive black holes 
(BHs) has been devoted to measuring their masses. This lead to the finding of the 
tight Mbh — f relation, which indicates the existence of strong links between the 
formation of the BH and of their host spheroids. Many scenarios have been proposed 
to explain this relation, and all agree on the key role of BHs' growth and feedback in 
shaping their host galaxies. However the currently available observational constraints, 
essentially BH masses and galaxy photometry, are not sufficient to conclusively select 
among the alternatives. A crucial piece of information of the black hole formation is 
recorded in the orbital distribution of the stars, which can only be extracted from high- 
resolution integral-field (IF) stellar kinematics. The introduction of IF spectrographs 
with adaptive optics on large telescopes opens a new era in the study of BHs by finally 
allowing this key element to be uncovered. This information will be complementary 
to what will be provided by the LISA gravitational wave satellite, which can directly 
detect coalescing BHs. Here an example is presented for the recovery of the orbital 
distribution in the center of the giant elliptical galaxy M87, which has a well resolved 
BH sphere of influence, using SAURDN IF kinematics. 

PACS numbers: 98.10.-I-7, 98.52.Eh, 98.62.Js 

Submitted to: Class. Quantum Grav. 22, (2005) S347-S353 
1. Introduction 

In recent years there has been tremendous progress, primarily from space-based 
observations, in our understanding of the distribution of central black hole (BH) masses 
and the relation between BHs and their host galaxies (see de Zeeuw 2004 for a review). 
The resulting picture of BH demography is summarized by the correlation between 
BH mass and absolute spheroid luminosity (e.g., Magorrian et al. 1998) and the much 
tighter correlation between BH mass and galaxy central velocity dispersion (the M^^ — a 
relation; Ferrarese & Merrit, 2000; Gebhardt et al. 2000). It is now clear that BHs are 
nearly ubiquitous in galaxies, and there is increasing evidence that BH formation plays a 
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Figure 1. Observable signatures of tiie formation of binary BHs in galaxies, predicted 
by iV-body simulations, (a) projected surface brightness. From top to bottom the 
different lines indicate the evolution, starting from the time th of the formation of 
the hard black hole binary. The vertical arrow indicates the location of the core 
radius Rc for the lowest profile, (b) orbital distribution (anisotropy; here defined as 
/3 = 1 — af/a"^, where cr,. is the radial and ut the tangential component of the velocity 
dispersion tensor) at the time th. The thick and thin lines are measured along the 
major and minor axis respectively. The dashed vertical line indicates the radius i?BH 
of the BH sphere of influence, while the solid vertical line corresponds to the radius 
Rc- (c) same as in the previous panel, but at the end of the simulation. Inside the 
radius where the density profile is 'scoured' by the BH binary, the orbital distribution 
becomes tangentially biased. For comparison the dotted line shows the anisotropy for 
a simulation with a single BH, where the orbital distribution remains mildly radially 
anisotropic. (Adapted from Milosavljevic & Merritt 2001). 



key role in driving galaxy formation, via feedback processes. The correlations obtained 
from nearby galaxies are also being extended to study the high-redshift universe. To 
date, however, the only observable constraints on the BH formation process are the BH 
masses and the nuclear stellar density profiles of nearby galaxies. Until now, very little 
is known about the nuclear orbital distribution of the stars, which constitutes a fossil 
record of the BH formation process. The advent of high spatial resolution integral-field 
spectroscopic observations is opening a new era in the study of BHs by allowing this 
crucial piece of information to be reliably extracted. 

2. Observable Constraints on Black Hole Formation Scenarios 

Spheroidal components of galaxies are thought to form through a series of hierarchical 
mergers of smaller systems, which themselves have central BHs. As two systems merge, 
their central BHs rapidly fall to the minimum of the potential well due to dynamical 
friction, and create a binary system at the center of the remnant. The binary looses 
energy and shrinks by capturing stars on radial orbits which pass nearby, and ejecting 
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them at much higher velocities, thus dynamically heating the central regions and 
lowering its stellar density, scouring out a 'core' of radius in the photometric profile 
(Quinlan & Hernquist 1997; Milosavljevic & Merritt 2001). In this way, a coalescing 
BH binary can eject a total mass of stars similar in order to its own mass, and influence 
the orbital structure of the remnant nucleus out to the radius Rc, which is several times 
larger than the actual BH radius of influence. The imprint of this process is found in the 
central orbital structure, which can be quantifled e.g. in terms of the ratio of radial to 
tangential orbits (Fig. 1). BHs can also grow by acquiring gas which sinks via dissipation 
to the bottom of the galaxy potential well. In this case, simulations show that a cusp 
forms in the density proflle when the gas subsequently forms stars, while the orbital 
distribution is only weakly affected. These two black hole formation mechanisms have 
been proposed as an explanation for the dichotomy of central luminosity proflles of early- 
type galaxies (Faber et al. 1997). While the details of the theoretical predictions for the 
effects of the BH growth on the orbital distribution are still under lively debate (Merritt 
& Milosavljevic 2004; Makino & Funato 2004), due to the difficulty of constructing 
accurate A'"-body simulations with very large numbers of stars near the BH singularity, 
it already appears clear that BHs have a profound infiuence on the structure of galaxy 
nuclei, which retain the fossil record of the formation process. One of the strongest 
observable constraints on the BH formation scenarios is its effect on the central orbital 
structure. 

3. Limitations of high-resolution long-slit studies 

Probing the BH radius of infiuence demands sub-arcsecond spatial resolution, even for 
nearby galaxies. For this reason, STIS onboard HST has been used for many well- 
determined BH masses measured from dynamical modeling of the kinematics of stars 
or gas. There are also several studies currently in progress to establish black hole 
masses using adaptive optics (AO)-assisted long-slit measurements from the ground, 
using NAOS-CONICA (VLT) and NIRI-Altair (Gemini). There are strong limitations, 
however, to what can be achieved by long-slit studies. A single slit observation is not 
sufficient to detect possible signs of non-axisymmetry in the kinematics, resulting from 
a nuclear bar, for example. Therefore the symmetry assumptions generally made in the 
dynamical models, used to measure BH masses, cannot be tested. 

More importantly, the orbital distribution (e.g., the anisotropy) cannot be recovered 
without knowledge of the line-of-sight velocity-distribution (LOSVD) of the stars at 
all spatial positions on the galaxy image projected on the sky (Fig. 2). Integral-field 
kinematic data is therefore the only way to reliably constrain the orbital structure 
(Cappellari et al. 2004), and general orbit-based methods can then be used in realistic 
conditions, for a given potential, to extract the orbital distribution from the two- 
dimensional (2D) kinematics (Krajnovic et al. 2005). Ignorance of the anisotropy due 
to only having a single long-sht directly translates into large uncertainties in the BH 
mass determination (Verolme et al. 2002), as well as invalidating any constraint on 




Figure 2. Each column from left to right shows the kinematic moments: mean 
velocity, a, and higher Gauss-Hermite moments h3-h6 (see text). Top row: 2D-binned, 
symmetrized and linearly interpolated SAURON observations of the early-type galaxy 
NGC4473 from Emsellem et al. (2004). Middle row: Three-integral Schwarzschild 
model fit using the full SAURON field. The analysis of this model shows that the galaxy 
contains two flattened countcrrotating spheroids, of unequal mass (Cappellari et al. 
2004). This explains the unusual high velocity dispersion along the major axis and the 
V-shaped velocity field. Bottom row: The same model, but fitted only to a simulated 
major-axis long slit extracted from the SAURON field. All models are regularized and 
have a smooth phase space. The kinematics are fitted perfectly along the slit, and 
the galaxy density is everywhere well reproduced; outside the slit, however, the model 
is completely different from the true galaxy. Only using the full 2D field are the 
kinematics, and therefore the orbital structure, properly constrained. 



the BH formation scenario. Although observational studies of the velocity anisotropy 
around BHs exist (Gebhardt et al. 2003), these are based on single-slit information, from 
HST/STIS, and so cannot place tight constraints on the orbital structure. Avoiding 
these issues by obtaining multiple slit positions has a prohibitive observational expense, 
and no such studies exist. 

Until now, no instrument could provide integral-field kinematics at the resolution 
of HST (but see Bacon et al. 2001a), so it was impossible to reliably infer the orbital 
distribution near the nuclear black holes. The advent of AO-fed IF units on large 
telescopes, like SINFONI (VLT), OASIS (WHT), NIFS (coming on Gemini) and OSIRIS 
(coming on the Keck telescopes), is finally going to change this situation. In particular, 
the future prospect of laser guide star-based AO systems on these instruments will allow 
statistically meaningful samples to be obtained. Observations from such instruments, 
combined with general dynamical models, will provide the first reliable determination 
of the effect of a black hole on the nuclear stellar orbits, giving a unique insight into the 
formation of BHs and their host systems. 
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Figure 3. Each column from left to right shows the kinematic moments: mean 
velocity, cr, and higher Gauss-Hermite moments hs-hg. Top row: Voronoi 2D- binned 
and symmetrized SAURDN observations of the giant elliptical galaxy M87. Bottom row: 
Three-integral Schwarzschild model fit to the above kinematics. 



4. Nuclear orbital distribution of M87 from integral-field kinematics 

To demonstrate what IF observations can provide, we present here the recovery of the 
orbital distribution from general axisjTiimetric three-integral dynamical modeling of the 
IF stellar kinematics of M87, which is part of the SAURON survey (Bacon et al. 2001b; 
de Zeeuw et al. 2002). This galaxy is an ideal candidate for this ground-based study 
as its BH sphere of influence Rbk ~ 1-4" is resolved even in natural seeing. Moreover, 
the core, whose structure is thought to have been shaped by the BH binary, has a 
radius Rc ~ 7.5" (Faber et al. 1997) and is very well sampled by the observed SAURON 
kinematics. 

The observations were presented in Emsellem et al. (2004) and we refer to that paper 
for details. In brief, the SAURON IF observations were spatially binned to a minimum 
S/N ^ 60 using the Voronoi 2D-binning algorithm by Cappellari & Copin (2003) and the 
stellar kinematics were subsequently extracted with the penalized pixel-fitting (pPXF) 
method of Cappellari & Emsellem (2004). This provided, for each Voronoi bin, the mean 
velocity V, the velocity dispersion a and the higher order Gauss-Hermite moments of 
the velocity, including h^-h^ (van der Marel & Franx 1993; Gerhard 1993). For the 
three-integral dynamical modeling, we adopted Schwarzschild's (1979) numerical orbit- 
superposition method, which can fit all kinematical and photometric observations (Rix 
et al. 1997; van der Marel et al. 1998; Cappellari et al. 2002). A similar approach was 
adopted by other groups to measure the black hole (BH) masses in galaxy nuclei (e.g. 
Gebhardt et al. 2003; Valluri et al. 2004). 

Our best fitting model with constant stellar mass-to-light ratio M/L (Fig. 3; for 
an assumed inchnation i = 90°) is able to reproduce all the details of the observed 
kinematics. As our kinematic measurements do not have enough spatial resolution to 
provide tight constraints on the BH mass, we assume Mbh ~ 3 x IO^Mq as derived from 
the HST gas kinematics (Harms et al. 1994; Macchetto et al. 1997). The anisotropy 
profile of the model is presented in Fig. 4. It shows that the main galaxy body is mildly 
radially anisotropic even inside Rc, and it is not until well inside Rbr (at the limit of 
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Figure 4. Radial profile of the orbital anisotropy Urjcrt in M87. The different 
solid lines are measured along different polar angles, in the galaxy meridional plane, 
and provide an indication of the model uncertainties. The dotted horizontal line 
corresponds to an isotropic model. The dashed vertical line indicates the radius i?BH 
of the BH sphere of influence, while the solid vertical line corresponds to the location 
of the break radius i?c, which is indicative of the core size. The main galaxy body is 
mildly radially anisotropic until well inside the sphere of influence of the BH (at the 
limit of our spatial resolution), where the orbits become tangentially anisotropic. 



our spatial resolution) that it becomes tangentially anisotropic. This result is consistent 
with previous models by Dressier & Richstone (1990), van der Marel (1994) and Merritt 
& Oh (1997) which only made use of the velocity dispersion along a single slit position. 
This is not surprising given the near circular symmetry of the projected observables 
of this galaxy in the central regions. In a nonrotating spherical stellar system, if the 
potential is given, the anisotropy is uniquely specified by the velocity dispersion alone 
(Binney & Manon 1982). Therefore, the fact that our model is able to fit the higher 
order moments of the velocity as well, gives us confidence that the orbital distribution 
is being correctly recovered. 

Although little can be learned from a single object and a specific simulation, it is still 
interesting to see how the observations and the predictions compare. We find that the 
observed distribution of M87, which is the prototypical giant core elliptical galaxy, is not 
obviously explained by current simulations of the BH binary core-scouring mechanism. 
In fact, although the observed anisotropy is consistent with the predicted one at the time 
th of the formation of the hard BH binary (Fig. lb), the steep surface brightness profile 
(Fig. la) is inconsistent with the observations (Fig. 7 of van der Marel et al. 1994). The 
observed surface brightness profile is better reproduced by the simulations at later times, 
when the binary had time to scour a deeper core, but the anisotropy profile at that time 
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(Fig. Ic) is inconsistent with the observed one. This discrepancy may be due to the 
effect of muhiple merger events (Milosavljevic & Merritt 2001), making the simulations 
and the observations difficult to compare in detail. This practical example demonstrates 
the method, and highlights the need to perform this kind of study for larger samples 
of objects, in order to tightly constrain the BH formation mechanism of these massive 
(Mbh> 10^) BHs in nearby galaxies, and uncover fossil signatures of BHs binaries, which 
can be compared with detailed predictions from simulations. This information will be 
complementary to what will be provided by the LISA gravitational wave satellite, which 
can directly detect coalescing BHs (if this happens) in the smaller mass range Mbh<10^ 
Mq, and in a much larger redshift range (up to 2;~10). Understanding the BHs growth is 
expected to shed light on the feedback processes by which BHs shape galaxy evolution. 
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